Coupling nonlinear Stokes and Darcy flow using mortar finite elements

We study a system composed of a nonlinear Stokes flow in one subdomain coupled with a nonlinear porous medium flow in another subdomain. Special attention is paid to the mathematical consequence of the shear-dependent fluid viscosity for the Stokes flow and the velocity-dependent effective viscosity for the Darcy flow. Motivated by the physical setting, we consider the case where only flow rates are specified on the inflow and outflow boundaries in both subdomains. We recast the coupled Stokes–Darcy system as a reduced matching problem on the interface using a mortar space approach. We prove a number of properties of the nonlinear interface operator associated with the reduced problem, which directly yield the existence, uniqueness and regularity of a variational solution to the system. We further propose and analyze a numerical algorithm based on mortar finite elements for the interface problem and conforming finite elements for the subdomain problems. Optimal a priori error estimates are established for the interface and subdomain problems, and a number of compatibility conditions for the finite element spaces used are discussed. Numerical simulations are presented to illustrate the algorithm and to compare two treatments of the defective boundary conditions.     

气体钻水平井环空冲蚀能量分析

在建立水平井流动方程的基础上,通过气固两相流动基本数学模型的求解数值模拟不同工况条件下的气固混合物冲蚀能量,得出了冲蚀能量的变化规律。算例数值模拟显示:环空速度剖面主要由注气量大小决定,机械钻速的影响不大;相同注气量条件下,机械钻速降低与井下岩屑浓度等比例下降;注气量的增加能提高井眼净化程度,但不如降低机械钻速效果明显;环空总冲蚀能中气体冲蚀能量占绝对主要地位,岩屑冲蚀能量只占极小部分,在相同注气量条件下总冲蚀能量基本相同;岩屑冲蚀能量主要取决于机械钻速,机械钻速越大井底产生岩屑越多冲蚀能力越强,相同机械钻速条件下,较大注气量能提供较大岩屑运移速度,提高了岩屑冲蚀能量。     

Radiometrie Investigations into Gas-Solid-Fluidized Beds

Investigations made into a gas-solid-fluidized bed model reactor, aiming at obtaining statements about the spatial and time density distributions within the fluidized bed and the correlation between the state of the turbulent layer and a characteristic measured quantity as well as variation of a physical property of the solid. statisticl evaluation techniques are indispensable, owing to the elementary actions occuring within the fluidized bed reactor and their random nature. Such parameters of the fluidized bed could be determined from a totality of individual measurements for measuring periods of 10 ms per single measurement by use of a miniature measuring probe within the interior of the fluidized bed reactor. This probe consisted of a beta-radition source and a Si-surface barrier-layer detector.     

A hybrid DSMC/Navier–Stokes frame to solve mixed rarefied/nonrarefied hypersonic flows over nano‐plate and micro‐cylinder

We extend a hybrid DSMC/Navier–Stokes (NS) approach to unify the DSMC and the NS simulators in one framework capable of solving the mixed non‐equilibrium and near‐equilibrium flow regions efficiently. Furthermore, we use a one‐way state‐based coupling (Dirichlet–Dirichlet boundary‐condition coupling) to transfer the required information from the continuum region to the rarefied one. The current hybrid DSMC–NS frame is applied to the hypersonic flows over nanoflat plate and microcylinder cases. The achieved solutions are compared with the pure DSMC and NS solutions. The results show that the current hybrid approach predicts the surface heat transfer rate and shear stress magnitudes very accurately. Some important conclusions can be drawn from this study. For example, although the shock wave region would be a non‐equilibrium region, it is not necessary to use a pure DSMC simulator to solve it entirely. This is important when the researchers wish to predict the surface properties such as velocity slip, temperature jump, wall heat flux rate, and friction drag magnitudes accurately. Our investigation showed that our hybrid solution time would be at least 40% (for the flat plate) and 35% (for the cylinder) of the time that must be spent by a pure DSMC solver to attain the same accuracy.Copyright © 2013 John Wiley & Sons, Ltd.     

Film thickness measurements in liquid–liquid slug flow regimes

At present there is significant interest in the development of small scale medical diagnostic equipment. These devices offer faster processing times and require smaller sample volumes than equivalent macro scale systems. Although significant attention has been focused upon their outputs, little attention has been devoted to the detailed fluid mechanics that govern the flow mechanisms within these devices. Conventionally, the samples in these small scale devices are segmented into distinct discrete droplets or slugs which are suspended in an organic carrier phase. Separating these slugs from the channel wall is a very thin film of the organic carrier phase.The magnitude of this film is the focus of the present study and the effects of sample slug length and carrier phase fluidic properties on the film are examined over a range of Capillary numbers. A non-intrusive optical technique was used to capture images of the flow from which the magnitude of the film was determined.The experimental results show that the film is not constant along the length of the slug; however above a threshold value for slug length, a region of constant film thickness exists. When compared with existing correlations in the literature, the experimental data showed reasonable agreement with the Bretherton model when the Capillary number was calculated based on the mean two phase flow velocity. However, significant differences were observed when the Capillary number was redefined to account for the mean velocity at the liquid interface, i.e., the mean slug velocity.Analysis of the experimental data revealed that it fell into two distinct flow regimes; a visco-capillary regime and a visco-inertial regime. A modified Taylor expression is presented to estimate the magnitude of the film for flows in the visco-capillary regime while a new model is put forward, based on Capillary and Weber numbers, for flows in the visco-inertial regime. Overall, this study provides some novel insights into parameters, such as aqueous slug length and carrier phase fluidic properties, that affect the thickness of the film in liquid–liquid slug flow regimes.     

FRACTURE ASSESSMENT OF AN INTERFACE CRACK BETWEEN TWO DISSIMILAR MAGNETOELECTROELASTIC MATERIALS UNDER HEAT FLOW AND MAGNETOELECTROMECHANICAL LOADINGS

A magnetoelectrically permeable interface crack between two semi-infinite magnetoelectroelastic planes under the action of a heat flow and remote magnetoelectromechanical loadings is considered, where the assumption of frictionless contact between two dissimilar halfplanes is adopted. Not only the solutions of the interface crack problem are presented in an explicit form, but also the general condition for the transition from a perfect thermal contact of two magnetoelectroelastic bodies to their separation is given.     

Surrogate-based modeling and dimension reduction techniques for multi-scale mechanics problems

Successful modeling and/or design of engineering systems often requires one to address the impact of multiple "design variables" on the prescribed outcome.There are often multiple,competing objectives based on which we assess the outcome of optimization.Since accurate,high fidelity models are typically time consuming and computationally expensive,comprehensive evaluations can be conducted only if an efficient framework is available.Furthermore,informed decisions of the model/hardware’s overall performance rely on an adequate understanding of the global,not local,sensitivity of the individual design variables on the objectives.The surrogate-based approach,which involves approximating the objectives as continuous functions of design variables from limited data,offers a rational framework to reduce the number of important input variables,i.e.,the dimension of a design or modeling space.In this paper,we review the fundamental issues that arise in surrogate-based analysis and optimization,highlighting concepts,methods,techniques,as well as modeling implications for mechanics problems.To aid the discussions of the issues involved,we summarize recent efforts in investigating cryogenic cavitating flows,active flow control based on dielectric barrier discharge concepts,and lithium(Li)-ion batteries.It is also stressed that many multi-scale mechanics problems can naturally benefit from the surrogate approach for "scale bridging."     

Cell face velocity alternatives in a structured colocated grid for the unsteady Navier–Stokes equations

The use of a colocated variable arrangement for the numerical solution of fluid flow is becoming more and more popular due to its coding simplicity. The inherent decoupling of the pressure and velocity fields in this arrangement can be handled via a special interpolation procedure for the calculation of the cell face velocity named pressure‐weighted interpolation method (PWIM) (AIAA J. 1983; 21 (11):1525–1532). In this paper a discussion on the alternatives to extend PWIM to unsteady flows is presented along with a very simple criterion to ascertain if a given interpolation practice will produce steady results that are relaxation dependent or time step dependent. Following this criterion it will be shown that some prior schemes presented as time step independent are actually not, although by using special interpolations can be readily adapted to be. A systematic way of deriving different cell face velocity expressions will be presented and new formulae free of Δt dependence will be derived. Several computational exercises will accompany the theoretical discussion to support our claims. Copyright © 2009 John Wiley & Sons, Ltd.     

Evaluation of linear and nonlinear convection schemes on multidimensional non‐orthogonal grids with applications to KVLCC2 tanker

Almost all evaluations of convection schemes reported in the literature are conducted using simple problems on uniform orthogonal grids; thus, having limited contribution when solving industrial computational fluid dynamics (CFD), where the grids are usually non‐orthogonal with distortions. Herein, several convection schemes are assessed in uniform and distorted non‐orthogonal grids with emphasis on industrial applications. Linear and nonlinear (TVD) convection schemes are assessed on analytical benchmarks in both uniform and distorted grids. To evaluate the performance of the schemes, four error metrics are used: dissipation, phase and L₁ errors, and the schemes' effective order of accuracy. Qualitative and quantitative deterioration of these error metrics as a function of the grid distortion metrics are investigated, and rigorous verifications are performed. Recommendations for effective use of the convection schemes based on the range of grid aspect ratio (AR), expansion ratio (ER) and skewness (Q) are included. A ship hydrodynamics case is studied, involving a Reynolds averaged Navier–Stokes simulation of a bare‐hull KVLCC2 tanker using linear and nonlinear convection schemes coupled with isotropic and anisotropic Reynolds‐stress (ARS) turbulence models using CFDShip‐Iowa v4. Predictions of local velocities and turbulent quantities from the midships to the nominal wake plane are compared with experimental fluid dynamics (EFD), and rigorous verification and validation analyses for integral forces and moments are performed for 0° and 12° drift angles. Best predictions are observed when coupling a second‐order TVD scheme with the anisotropic turbulence model. Further improvements are observed in terms of prediction of the vortical structures for 30° drift when using TVD2S‐ARS coupled with DES. Copyright © 2009 John Wiley & Sons, Ltd.